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We report a search for B decays to selected final states with the rjc me¬ 
son: B^ —>• , B^ —)■ K^rjcUJ, B^ —> and B^ —)■ 

The analysis is based on 772 x 10® BB pairs collected at the T(4S') reso¬ 
nance with the Belle detector at the KEKB asymmetric-energy e+e" col¬ 
lider. We set 90% conhdence level upper limits on the branching fractions 
of the studied B decay modes, independent of intermediate resonances, in 
the range (0.6 —5.3) x 10“^. We also search for molecular-state candidates 
in the D^D^ + D^D^ and D*^D*\d*°D*^ combinations, 

neutral partners of the Z(3900)^ and Z(4020)^, and a poorly understood 
state A(3915) as possible intermediate states in the decay chain, and set 
90% conhdence level upper limits on the product of branching fractions to 
the mentioned intermediate states and decay branching fractions of these 
states in the range (0.6 — 6.9) x 10“®. 
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1 Introduction 


The state X(3872) was first observed by Belle in exclusive J/'ip de¬ 

cays p. The X(3872) mass is close to the -|- m^*o threshold, which engen¬ 
dered a hypothesis that this state may be a D^D*^ molecule [2]. Angular analysis of 
the X(3872) —> decay by LHCb |3j determined all its quantum numbers: 

= 1++. If X(3872) is indeed a D^D*^ molecule, there can exist other “X(3872)- 
like” molecular states with different quantum numbers. Some may reveal themselves 
in the decays to final states containing the r]c meson. For example, a D^D*^ — D^D*^ 
combination (denoted hereinafter by Xi(3872)) with quantum numbers 
would have a mass around 3.872 GeV/c^ and would decay to rj^p and Combi¬ 
nations of D^lf +lfD^, denoted by X(3730), and D*^D*^ + D*^D*^, denoted by 
X(4014), with quantum numbers = O’*"’*' would decay to rjcV and . The 
mass of the X(3730) state would be around 2mD0 = 3.730 GeV/c^ while that of the 
X(4014) state would be near 2mD*o = 4.014 GeV/c^. 

Recently, a new charged state Z(3900)^ was found in y(4260) decays by Belle [1] 
and BESIII p. Since this particle is observed in the decay to it should 

contain at least four quarks. The Z(3900)^ was confirmed in the decay to 
by an analysis of CLEO-c data [6] that also reported evidence for its neutral isotopic 
partner ^(3900)*^. Another exotic charged state Z(4020)^ was observed by BESIII in 
decays to P and {D*D [8]. There are some indications from these analyses 
that the spin and parity of the charged states might be = I"*". 

The near-threshold enhancement in the uJ/ij) invariant mass distribution named 
X(3940) was first observed by Belle in exclusive B —)■ KojJ/'ip decays [9]. Later, in 
the same decay mode, BaBar discovered X(3915) [10]. The parameters of y(3940) 
are consistent with those of X(3915), so they are considered to be the same particle. 
The quantum numbers of X(3915) are claimed to be = O’*"'*', but its nature is 
still undetermined, and there are several interpretations describing this state [in- 

To search for the particles described above, we reconstruct 'Pc mesons via the 
KgK^'ir'^ mode and study the following four decays of charged B mesons: 

1. the (tt+tt”) decay mode: B^ —>■ K^X —>■ where we look for 

Xi(3872), Z(3900)° and Z(4020)°; 

2. the (cu) decay mode: B^ —)■ K^X —)■ where we look for Xi(3872); 

3. the {rj) decay mode: B^ —>■ K^X —)■ K^^rjcV), where we look for X(3730), 
X(4014) and X(3915); 

4. the (vr®) decay mode: B^ —>• K^X —>■ X^(uc7r°), where we look for X(3730), 
X(4014) and X(3915). 

The analysis is based on a data sample that contains 772 x 10® BB pairs, collected 
with the Belle detector at the KEKB asymmetric-energy e’''e“ collider. 
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2 ^ K'^Vc + hadrons 

We search for D(*)0£)(*)o molecular-state candidates Z(3900)°, Z(4020)°, and X(3915) 
in the following B meson decays: —)■ K'^r]cTT~^7r~, B^ K'^rjcUJ, —)■ 

and —)■ K^rjcTT^. To determine the branching fractions, we perform a binned 

maximum-likelihood fit of the AE distribution that is modelled by a peaking signal 
and featureless background. For the ( 7 r’'' 7 r“), (cn) and (tt^) decay modes, the signal 
function is the sum of two Gaussians {G) and the background function is a linear 
polynomial. Here and in the following the detector resolution is taken into account. 
The AE distribution for the (cn) mode is shown in Figured] (left). In the (tt+tt”) and 
(tt^) modes, we observe some significant signal and so perform a two-dimensional fit 
of the KsKtt invariant mass and AE distributions shown in Figure [2j Since the rj 
candidate is reconstructed in two decay modes, we perform a combined fit of the AE 
distribution corresponding to 77 —)■ 77 and 7 —)■ 7i~^7r~7T^ shown in Figure [1] (middle 
and right). The fit results are summarized in Tabled! 





Figure 1: The AE distribution for the decay —>■ K^rjcOJ (left). The combined fit 
projections of the AE distributions in case of the 7 —>■ 77 (middle) and 7 —>■ Tr+vr^Tr^ 
(right) modes for the decay B^ —>■ 


Decay mode 

Efficiency, % 

Yield 

—> K^rjcUJ 

0.53 ±0.01 

-41 ± 27 

B^ —> K^r]cT^^Ti~ 

2.84 ±0.02 

155 ± 72 

B^ K^r]c7r^ 

—> K^rjcT], 

3.69 ±0.01 

-1.9 ± 12.1 

7 —>• 77 

3.05 ±0.01 

-14 ± 26 

Tj —>■ 7r'''7r“7r° 

0.69 ±0.01 

-1.8 ±3.4 


Table 1: Fit results for B decays independent of intermediate resonances. 
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Figure 2: Projections of the two-dimensional fit in KsKir invariant mass (left) and 
AE (right) for the decays —)■ K^t]cT^^'k~ (top) and —>■ (bottom). Each 

projection is plotted for events in the whole fitting range of the other projection. 

3 Xi(3872), JV(3730) and X(4014) 

The invariant mass distributions for (tt+tt”), (cu), [t]), and (vr^) modes are shown in 
Figures [3HS1 The corresponding signal yields are presented in Table [2J We validate 
our fit procedure by applying it to the decay K^'ip{2S), 'ip{2S) —)■ . 

This decay is similar to the (tt+tt”) decay except that we reconstruct the iIj{2S) 
meson in place of the Xi(3872) and the J/'ip m place of the 7]c- The J/ip meson, like 
the rjc, is reconstructed via the KgK^ir^ hnal state. We £t the J/'ip7r~^7T~ invariant 
mass distribution and obtain the number of signal events W = 20.2 ± 6.5, which 
corresponds to a signihcance of 3.5 standard deviations (cr). The expected number 
of events estimated using the world averages of the known branching fractions [12] is 
22 ± 4, which is consistent with Ng. 

4 Z(3900)° and ^(4020)° 

We perform a sequence of binned maximum likelihood fits of the r]c'n'~^7T~ invariant 
mass using the convolution of a Breit-Wigner and a Gaussian for the signal and a 
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Decay mode 

Efficiency, % Yield 

Xi(3872) —>• r^cVr+TT" 

7.95 ±0.02 17.9 ±16.5 

X]^(3872) —y 

1.92 ±0.02 6.0 ±12.5 

X(3730) ^ T],r], 


h ^ 77 

6.57 ±0.02 13.8 ±9.9 

r] —)■ 7r+7r“7r° 

1.18 ±0.01 1.4 ±1.0 

X(3730) ^ r]^TT^ 

6.52 ±0.02 -25.6 ±10.4 

X(4014) ^ 


7 ^ 77 

7.09 ±0.02 8.9 ±11.0 

Tj — )■ 7r^n~7r^ 

1.78 ±0.01 1.3 ±1.6 

X(4014) ^ 

7.55 ±0.02 -8.1 ±13.2 

X(3915) T]cT], 


T] ^ 'f'f 

6.60 ±0.02 -7.4 ±14.5 

Tj — )■ 7r’''7r“7r° 

1.64 ±0.01 -1.1 ±2.1 

X(3915) —>• TjcTT^ 

6.88 ±0.02 -4.3 ±18.1 


Table 2: Fit results for the Xi(3872), X(3730), X(4014) and X(3915) resonances. 


linear polynomial for the background. The Gaussian models the detector resolution, 
which is assumed to be similar to that obtained in Ref. na and equal to 9.8 MeV/c^. 
The Breit-Wigner mass is conhned to a 20 MeV/c^ window (the so-called mass bin) 
that is scanned in 20 MeV/c^ steps across the range (3.79 — 4.01) GeV/c^ for the 
Z(3900)° and (3.93 — 4.07) GeV/c^ for the Z(4020)°. The width is fixed to the 
weighted mean of the previously measured values (35 MeV/c^ for the Z(3900)° and 
12 MeV/c^ for the Z(4020)°). The obtained signal yield is shown in Figure El 

5 X(3915) 

The rjcT] and tJcTT^ invariant mass distributions are shown in Figures [7] and El The £t 
results are summarized in Table [2l 
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3.7 3.8 3.9 4 

M(ti,co), GeV/c2 


Figure 3: The distributions of r^cvr+vr (left) and t]cUJ (right) invariant masses in the 
search for the Xi(3872). 






Figure 4: The combined fit projections of the rjcr] invariant mass distributions in case 
of the 7] ^ 77 (left) and rj —)■ (right) modes corresponding to the search for 

the X(3730) (top) and X(4014) (bottom) resonances. 
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Figure 5: The r/cVr® invariant mass distributions corresponding to the search for the 
X(3730) (left), X(4014) (middle) andX(3915) (right) resonances. 



M(T|r7t+7r-), GeV/c^ 



M(T|r7t+7r-), GeV/c^ 


Figure 6: Dependence of the signal yield of the Z(3900)° (left) and Z(4020)° (right) 
on the mass bin. The mass bin is a 20 MeV/c^ window to which the mass is confined 
and scanned in 20 MeV/c^ steps across the fit range. 




Figure 7: The combined fit projections of the invariant mass distributions in case 
of the ?7 ^ 77 (left) and r] —> vr+Tr^Tr® (right) modes corresponding to the search for 
the X(3915) resonance. 
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6 Results 


Our data sample does not permit us to measure the branching products of production 
and decay of the states listed above nor the B decay branching fractions so we set 
upper limits instead, taking into account the statistical and systematic uncertainties. 
Upper limits on the branching fractions and products for all the studied decay modes 
are shown in Tables [3] and HI 


Decay mode 

Upper limit (90% G.L.) 


3.9 X 10"^ 

B^ K^7]cU: 

5.3 X 10-^ 

B^ —>• K^TJcTj, 

2.2 X 10“^ 

R± ^ K^r]c7r^ 

6.2 X lO-*^ 


Table 3: Results of branching fraction measurements for the B decays without an 
intermediate resonance. 


Resonance Decay mode 

Upper limit (90% G.L.) 

Xi(3872) r]c7r^n~ 

3.0 X 10"^ 

r]cUJ 

6.9 X 10“^ 

X(3730) peh 

4.6 X 10"^ 


5.7 X 10-® 

X(4014) rjcT] 

3.9 X 10"^ 


1.2 X 10“^ 

Z(3900)° r^cVr+TT” 

4.7 X 10-^ 

Z(4020)° 

1.6 X 10“^ 

X(3915) 

3.3 X 10-^ 

r]c7T^ 

1.8 X 10“^ 


Table 4: Results of branching fraction measurements for the B decays containing 
an intermediate exotic resonance. For Z(3900)° and Z(4020)° resonances the results 
are shown under the assumption that the masses are close to those of their charged 
partners. 


In case of the Z(3900)° and Z(4020)° mass scan, no signihcant signal is seen in 
any of the invariant mass bins. For the Z(3900)° resonance, we set upper limits 
in the range (1.8 — 4.7) x 10“® for the mass region (3.79 — 4.01) GeV/c^. For the 
^(4020)*^ resonance, we set upper limits in the range (1.6 — 3.7) x 10“® for the mass 
region (3.93 — 4.07) GeV/c^. If we assume that the Z(3900)° and Z(4020)° masses 
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are close to those of their charged partners, we obtain the upper limits on the product 
branching fractions shown in Table 01 
The results are published in [T4] . 
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